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Abstract 

The  effect  of  a  thin  viscous  fluid-mud  layer  on  nearshore  nonlinear  wave-wave  interactions  is  studied  using  a  parabolic  frequency-domain 
nonlinear  wave  model,  modified  to  incorporate  a  bottom  dissipation  mechanism  based  on  a  viscous  boundary  layer  approach.  The  boundary-layer 
formulation  allows  for  explicit  calculation  of  the  mud-induced  wave  damping  rate.  The  model  performed  well  in  tests  based  on  laboratory 
data.  Numerical  tests  show  that  damping  of  high  frequency  waves  occurs,  mediated  by  “difference”  nonlinear  interactions.  Simulations  of 
2-dimensional  wave  propagation  over  a  mud  “patch”  of  finite  extent  show  that  the  wave  dissipation  causes  significant  down  wave  diffraction 
effects. 

©  2007  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  propagation  of  ocean  surface  waves  has  been  studied 
intensively  in  the  last  few  decades,  particularly  as  numerical 
models  of  wave  propagation  over  varying  bathymetry  have 
become  sufficiently  skillful  for  use  in  operational  forecasting, 
engineering  studies,  etc.  Both  phase-resolving  models  via  the 
mild-slope  (e.g.  Berkhoff,  1972)  or  Boussinesq  (e.g.  Peregrine, 
1967)  equations,  and  phase-averaged  models  (e.g.  Booij  ct  al., 
1 999),  have  been  extensively  validated  for  cases  in  which  the 
bottom  composition  is  not  a  significant  concern  to  the 
propagation  characteristics  (sandy  beaches).  Wave-bottom 
frictional  effects  are  typically  of  secondary  importance  in 
sandy  environments.  On  spatial  scales  typical  of  surf  zones  on 
sandy  beaches,  wave  breaking  tends  to  overwhelm  bottom 
friction  dissipation  rates.  Over  larger  scales  (100-km  wide 
North  Carolina  shelf),  swell  energy  loss  rates  of  50-75%  have 
been  observed,  attributed  to  the  interaction  of  waves  with 
bottom  ripples  (Ardhuin  et  al.,  2002;  Ardhuin  ct  al.,  2003). 

By  comparison,  the  strength  of  dissipation  processes  specific 
to  wave  propagation  in  fine-grained,  cohesive  sedimentary 
environments  is  well  known.  Extreme  dissipation  rates  have 
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been  routinely  reported:  80%  wave  energy  loss  over  only  2.6 
wavelengths  (Gade,  1958,  laboratory  experiments);  more  than 
90%  energy  loss  across  the  20-km  wide  shallow  mudflats  (Wells 
and  Coleman,  1981,  off  the  coast  of  Surinam);  95%  incident 
energy  loss  across  1.1  km-wide  mudbanks  (Mathew,  1992,  off 
the  coast  of  India),  and  so  on. 

Theoretical  formulations  of  bed-induced  wave  dissipation 
are  based  on  the  assumption  that  wave  motion  reaches  the 
bottom  and  interacts  directly  with  bed  sediments.  The 
approaches  typically  focus  on  a  single  (dominant)  physical 
mechanism,  defined  by  the  bed  state  and  composition.  A 
number  of  physical  mechanisms  for  wave  dissipation  over 
muddy  seabeds  have  been  proposed  over  time,  based  on 
different  models  of  sediment  rheology.  The  diversity  of  mud 
states  and  corresponding  beds  include  poro-elastic  solids 
(Yamamoto  et  al.,  1978;  Yamamoto  and  Takahashi,  1985), 
viscous  Newtonian  fluids  (Dalrymple  and  Liu,  1978;  Ng,  2000), 
Bingham  fluids  (Mei  and  Liu,  1987),  generalized  Voight  solids 
(Macpherson,  1980;  Hsiao  and  Shcmdin,  1980;  Jiang  and 
Mehta,  1995,  1996),  and  non-Newtonian  fluids  (Chou  ct  al., 
1993;  Foda  et  al.,  1993).  With  the  exception  of  fluidization 
processes  (Foda  et  al.,  1993;  DeWit,  1995),  these  models  focus 
on  a  single,  well-defined  mud  phase.  Although  the  models 
predict  greatly  differing  damping  effects  (Lee,  1995),  each  may 
be  applicable  in  some  situations  (Mei  and  Liu,  1987). 
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All  the  above  mechanisms  require  direct  interaction  between 
the  wave  motion  and  sediment  bed,  and  thus  are  essentially 
long-wave  dissipation  mechanisms.  On  time/spatial  scales 
comparable  to  those  of  wave  evolution  on  a  sandy  beach 
(102r  or  102A,  with  T  and  A  characteristic  wave  length  and 
period),  direct  wave-bed  interaction  is  arguably  the  dominant 
dissipation  mechanism  outside  the  surf  zone.  However,  many 
important  cohesive  sediment  coasts,  (e.g.  Louisiana  and 
Mississippi  Bights,  USA),  exhibit  morphological  characteristics 
which  allow  wave  evolution  and  dissipation  on  much  larger 
scales.  For  example,  on  the  inner  shelf  fronting  Atchafalaya 
(Louisiana,  USA),  the  10-m  isobath  is  in  some  places  approxi¬ 
mately  50  km  offshore.  Mediated,  indirect  wave-sediment 
energy  exchange  processes,  typically  also  second  order  in 
magnitude,  can  become  important  over  such  scales.  Examples 
are  short-wave  dissipation,  and  processes  at  the  interface 
between  a  fluid-mud  layer  and  the  upper  column,  such  as 
surface-interface  wave  interactions  (Foda,  1989;  Hill  and  Foda, 
1998;  Jamali  et  al.,  2003),  or  Kelvin-Helmholtz/Holmboe-type 
lutocline  instabilities  (e.g.  Mehta  and  Srinivas,  1993). 

A  series  of  recent  field  experiments  (Sheremet  and  Stone, 
2003;  Sheremet  et  al.,  2005)  on  the  Atchafalaya  shelf  have 
raised  questions  about  the  possible  role  of  surface  wave-wave 
interactions  in  mediating  short  wave  dissipation  in  cohesive 
sedimentary  environments.  The  observations  showed  unexpect¬ 
ed  short-wave  damping  in  areas  with  cohesive  bed  sediments,  in 
contrast  to  normal  (negligible)  short-wave  damping  over  areas 
with  sandy  bottoms.  Sediment  monitoring  devices  subsequently 
deployed  at  the  same  muddy  site  showed  that  strong  spectrum¬ 
wide  wave  damping  coincided  with  formation  of  a  fluid-mud 
layer  with  sediment  concentrations  of  over  10  kg/m2. 

Nonlinear  wave-wave  interaction  processes  should  be 
important  over  a  20-km  shallow  (15-20  m  depth)  muddy  shelf. 
On  a  typical  sandy  beach  (slopes  1%)  near-resonant  3-wave 
interactions  are  known  to  cause  significant  energy  transfers 
over  a  few  wavelengths  (Freilich  and  Guza,  1984;  Elgar  and 
Guza,  1985;  Agnon  et  al.,  1993;  Kaihatu  and  Kirby,  1995;  Agnon 
and  Sheremet,  1997,  2000;  Herbcrs  and  Burton,  1997  and  many 
others).  Strong  phase  correlations  build  up  between  spectral 
components,  and  the  wave  field  is  no  longer  Gaussian.  However, 
the  role  of  wave  nonlinearity  in  the  evolution  of  wave  spectra 
in  environments  characterized  by  strong  (frequency  dependent) 
damping  is  not  well  understood.  Jiang  and  Zhao  ( 1 989)  and  Jiang 
et  al.  (1990)  treated  the  problem  of  solitary  and  cnoidal  wave 
propagation  (respectively)  over  fluid-mud  seabeds,  using  several 
mud  layers  of  different  viscosities  to  enact  the  damping  and 
seeking  analytical  closed-form  solutions.  Sheremet  et  al.  (2005) 
incorporated  the  mud  dissipation  model  of  Jiang  and  Mehta 
( 1 995 )  into  the  nonlinear  1 D  triad  shoaling  model  of  Agnon  et  al. 
(1993)  and  tested  the  hypothesis  of  Sheremet  and  Stone  (2003) 
concerning  the  dampening  of  high  frequency  waves.  Their 
numerical  simulations  based  on  nonlinear  wave  propagation 
models  suggest  that  nonlinear  energy  transfer  within  the  wave 
spectrum  may  be  important,  perhaps  providing  the  coupling 
between  the  short-  and  long-wave  spectral  bands,  allowing  energy 
to  flow  toward  long  waves,  where  it  can  be  efficiently  dissipated 
via  direct  wave-bottom  interaction. 


In  this  study  we  incorporate  a  mud-induced  energy 
dampening  mechanism  into  a  phase-resolving  nonlinear  (one- 
dimensional  and  parabolic  two-dimensional)  frequency  domain 
model  for  surface  waves.  One  of  the  attractions  of  the  phase- 
resolving  approach  is  the  direct  correspondence  between  pro¬ 
perties  of  the  model  realizations  and  data  from  various  remote 
sensing  platforms  (e.g.  interferometric  synthetic  aperture  radar, 
surface-mapping  lidar,  video).  This  offers  great  potential  for 
data  assimilation  and  inverse  modeling  to  deduce  aspects  of  the 
bottom  topography  and  sediment  heterogeneity  over  a  large 
area.  We  first  compare  the  resulting  model  to  data  from  a  lab¬ 
oratory  experiment,  then  demonstrate  some  interesting  proper¬ 
ties  of  the  model  in  one-  and  two-dimensional  scenarios. 

The  model  developed  here  describes  mud-induced  surface 
wave  dissipation  and  nonlinear  processes  but  does  not  account 
for  the  complexities  of  the  coupling  between  sediment  and  wave 
states  (feedback  between  waves  and  fluid  mud  layer).  However, 
our  focus  at  the  present  is  on  a  wave  model  which  can  be 
applied  in  general  situations  of  wave  propagation  over  muddy 
bottoms,  rather  than  a  detailed  investigation  of  the  processes 
at  the  interface.  This  is  not  unlike  many  previous  studies  de¬ 
tailing  the  representation  of  wave  breaking  dissipation  into 
a  surface  wave  model  by  the  inclusion  of  an  additional  term 
(e.g.  Kirby  and  Dalrymple,  1986). 

2.  Viscous  dissipation  mechanism  of  Ng  (2000) 

The  dispersion  relation  for  a  linear  system  can  be  generically 
written  as: 

Q(a>,K,h)  =  0 ,  (1) 

which  relates  wavenumber  K  to  water  depth  h  and  frequency  co 
(Whitham,  1974).  For  linear  dissipative  systems  the  solutions  K 
of  the  dispersion  relation  (1)  are  complex,  and  are  typically 
found  by  a  laborious  search  in  the  complex  plane  (e.g.  the 
complex  secant  method  used  by  Dalrymple  and  Liu,  1978).  The 
accuracy  of  the  calculation  of  the  complex  roots  is  also  an 
important  concern,  particularly  when  incorporating  the  algo¬ 
rithms  into  a  numerical  model  with  nonlinear  terms  which  arc 
sensitive  to  the  nature  and  degree  of  dissipation,  as  has  been 
observed  with  surf  zone  breaking  (Kirby  and  Kaihatu,  1996). 
Numerical  applications  of  wave  modeling  require  an  efficient 
and  robust  approach  for  calculating  these  roots,  particularly  in 
the  case  of  nonlinear  models  which  spend  significant  compu¬ 
tational  resources  on  the  calculation  of  nonlinear  terms.  Here, 
we  investigate  the  simplified  approach  developed  by  Ng  (2000), 
which  has  the  potential  to  approximate  the  physics  of  mud- 
induced  dissipation  well  enough  to  allow  reliable  calculation  of 
the  dissipative  modes  K. 

Laboratory  observations  by  Gade  (1958)  and  numerical 
simulations  based  on  the  model  developed  by  Dalrymple  and 
Liu  (1978)  suggest  that  mud-induced  wave  dissipation  is 
maximized  when  the  mud  layer  is  approximately  1.5  times  the 
boundary  layer  depth  in  the  viscous  fluid,  thus  setting  up  the 
requisite  scales  for  a  boundary  layer  solution.  The  idea  was 
pursued  by  Ng  (2000),  who  formulated  the  viscous  dissipation 
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mechanism  as  a  boundary  layer  reduction  of  the  Dalrymple  and 
Liu  two-layer  fluid  model.  Assuming  that  the  depth  of  the  mud 
layer  is  comparable  to  the  mud  layer’s  Stokes  boundary  layer 
thickness  implies  that 

ka~kdm*kSm<K  1,  (2) 


The  solutions  ((4)— (5))  is  a  function  of  nondimensional  pa¬ 
rameters:  y -pjpm,  where  p  is  the  density;  ( = 6J6W = (vjvw) 1  a ; 
and  d=d/6m.  The  subscript  w  refers  to  water.  The  real  part  of  the 
second  order  correction  k2  represents  the  modification  to  the 
wavelength  due  to  the  viscous  mud.  The  total  real  part  of  the 
wavenumber  k  is  determined  as: 


where  k=Re(K ),  a  is  the  wave  amplitude,  dm  the  depth  of  the 
mud  layer,  and  S  =  yjlv/b  is  the  Stokes  boundary  layer  thick¬ 
ness,  where  v  is  the  kinematic  viscosity  and  a  is  the  angular 
frequency  of  the  wave.  The  subscript  m  denotes  mud.  As  a 
result  of  this  scaling,  the  wavenumber  of  a  dissipative  mode  can 
be  written  as 


K  =  k\  +  k2,  with  |*i|»|*2|» 


(3) 


where  the  leading  order  term  k\  is  real  and  satisfies  the  non- 
dissipative  dispersion  relation 


co2  =  g/r|tanh^|/i, 


(4) 


written  in  the  absence  of  any  fluid-mud  layer.  The  second  order 
correction  k2  is  complex  and  incorporates  the  effect  of  the  fluid 
mud.  For  k2 ,  Ng  (2000)  obtains 


*2  =  “ 


Bk\ 


(5) 


sinh  k\  h  cosh  k\  h  4-  k\  h 
The  imaginary  part  of  the  wavenumber  k2  is  the  damping  rate: 
Im(B)k\ 


__  \  _ _ 

sinh  k\h  cosh  k\hk\h 

where  B  is  a  complex  coefficient 


Im(S)  = 
with 


2*3 

k\dm(B\  +  B2) 
IBy  1 


(6) 

(7) 

(8) 


B\  =  y(— 2y2  +  2y  —  1  —  (2)sinht7  coshJ  -y2( 
x  (cosh 2  A  -fsinh2  A)  —  (-y2  —  l)2C 
x  (cosh2  A  cos2  A  -l-sinh2  A  sin2  A)  —  2y{  1  —  y) 
x  ((  cosh  d  +y  sinh  d )  cos  d  , 

B2  —  y{-2y2  -f  2y  -  1  4-  (2)sin  A  cos  A  -2y(l  -  y) 
x  ((sinhJ+y  cosh  A)  sin  A, 


(9) 


(10) 


By  =  ((cosh  A  -fysinh  A)1  cos2  A  +  ((sinh  A  -Fycosh^^sin2  A 


(11) 


k  —  k\  +  Re(k2)  =  k\ 


Re(B)k  i 

sinhiti  h  cosh^i  h  4-  k\  h 


(12) 


The  solution  ((4)— (5))  is  remarkable  in  that  mud  effects  are 
explicitly  expressed  in  terms  of  the  first  order,  nondissipative, 
real,  wavenumber  k\. 

The  characteristics  of  the  dissipation  (6)  had  been  dis¬ 
cussed  by  Ng  (2000,  see  Figs.  2  and  3);  however  the  effect  of 
a  thin  viscous  layer  of  mud  on  the  wavenumber  k  was  not 
similarly  explored  by  Ng  (2000).  This  companion  effect  of 
mud  on  the  real  part  k  of  the  wavenumber  is  important  for 
nonlinear  effects  in  wave  spectra  since  the  wavenumbers 
determine  the  amount  of  detuning  away  from  resonance,  and 
thus  directly  affect  energy  transfer  within  the  wave  spectrum. 
Fig.  1  shows  the  effect  of  the  mud  on  real  part  of  the 
wavenumber  for  two  different  values  of  C  (C=10)  and 
C=100),  where  the  density  ratio  y  =  0.9;  the  boundary  layer 
scale  6m  was  kept  constant  while  the  mud  layer  depth  dm  was 
varied  with  the  dimensionless  depth  A.  Clearly  the  increasing 
mud  layer  causes  a  decrease  in  the  wavenumber  k  relative  to 
that  with  no  mud.  The  convention  typically  used  in  theories  of 
wave  propagation  over  mud  has  the  mud  layer  adding  to  the 
total  water  depth;  the  wave  experiences  an  increase  in  water 
depth  with  the  addition  of  the  mud  layer,  and  the  resulting 
decrease  in  wavenumber  is  similar  to  an  inverse  shoaling 
effect.  The  degree  to  which  the  wavenumber  decreases  is 
dependent  on  the  mud  properties.  A  comparison  to  a  similar 
figure  shown  in  Dalrymple  and  Liu  (1978;  Fig.  3)  reveals  that 
the  present  theory  has  a  smaller  rate  of  wavenumber  decrease 
to  increase  of  d.  This  is  sensible  since  the  validity  of  the 
boundary  layer  approximation  in  Ng  (2000)  decreases  with 
increasing  d. 

Fig.  2  shows  the  variation  of  (6),  normalized  by  the 
maximum  dissipation  Dmax  for  each  parameter  set,  with  kh. 
This  normalization  is  done  to  help  visualize  the  shape  of  the 
frequency  dependence  and  reveal  the  similarity  across  mud 
parameter  sets;  the  value  of  Dmax  for  these  parameters  actually 
varies  over  three  orders  of  magnitude.  For  this  analysis,  the 
water  depth  h  =  \  m.  We  use  y  =  (0.5,0.9)  and  C=(10,100).  For 
c=  10  we  set  dm= 0.02  m,  and  for  (=100  we  use  dm-02  m. 
This  allows  d  to  have  the  same  values  over  the  range  of  kh  for 
all  mud  parameters  used.  Despite  the  differences  in  mud 
parameters,  the  curves  have  a  very  similar  shape  and  maxima 
which  are  located  near  kh  ~0.88.  All  curves  are  characterized 
by  a  rapid  rise  from  minimal  dissipation  for  kh  <0.88  and  a 
relatively  slow  fall  to  zero  damping  with  kh  >0.88.  This  is 
because  waves  in  very  shallow  water  (low  kh)  will  have 
mud  boundary  layer  thicknesses  Sm  which  are  greater  than 
the  depth  of  the  mud  layer  (low  A).  Relatively  little  damping 
thus  occurs.  Conversely,  waves  in  deeper  water  (high  kh)  will 
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Fig.  1.  Effect  of  viscous  mud  layer  on  wavenumber.  Solid  line:  C=  100,  7=0.9.  Dashed  line:  (  =  10,  7=0.9. 


not  be  as  influenced  by  the  bottom,  and  thus  also  experience 
little  damping.  The  damping  appears  to  be  at  least  33%  of 
the  maximum  value  in  the  range  0.1  <  kh  <2.4,  encompass¬ 
ing  the  upper  end  of  the  shallow  water  range  (kh<  tt/10) 
and  most  of  the  intermediate  depth  range  (tt/10  <  kh  <tt). 
We  note  here  that  different  values  of  dm  will  lead  to  a  differ¬ 
ent  range  of  d  over  the  range  of  kh  considered.  This  will  in 
turn  lead  to  different  values  of  kh  at  which  the  dissipation  is 
maximum. 

3.  Numerical  model  for  nonlinear  wave  transformation 
over  a  viscous  mud  layer 

We  now  turn  our  attention  toward  the  implementation  of  this 
mechanism  into  a  wave  model.  This  is  done  as  a  superposition 
of  an  added  dissipation  effect  into  the  model,  as  previously 
noted. 

We  use  the  nonlinear  model  of  Kaihatu  and  Kirby  (1995), 
which  is  a  parabolic  frequency-domain  mild-slope  equation 
model  with  second-order  nonlinear  wave-wave  interactions. 
The  free  surface  rj  is  assumed  to  be  expressible  as: 

ri(x,y,t)  =  ~  +  c.c  (13) 

where  An  is  the  complex  amplitude  of  the  free  surface  elevation, 
“c.c.”  denotes  complex  conjugate  of  the  term  in  the  square 
brackets,  and  kn  is  the  wavenumber  associated  with  the  radian 
frequency  con.  With  this  and  the  boundary  value  problem  for 
water  waves  (extended  to  second  order  in  wave  amplitude  a ), 


Kaihatu  and  Kirby  (1995)  derived  the  following  parabolic 
equation  for  An: 


dA 


2i{kCCg)-±  -  2(kCCg)n(k„  -kn)An  +  i 


d(kCCg) 


dx 


An 


k„)dx 

(14) 


where  C  is  the  phase  velocity,  Cg  is  the  group  velocity,  and  R 
and  5  are  interaction  coefficients 

RnJ  =  — - —  [(o2nk/kn-i  +  ( k,  +  k„-, ){a)n-, k,  +  c o, *„_,)«„] 

0)l(On-l 

O}2  ,  ,  ~i  \ 

- 1  ((Of  +  (O/OJn-l  +  (02_, )  (15) 

8 

Snj  =  +  {kn+i  —  ki)((on+iki  +  wikn+i)u)„] 

(OiU)n+i 

CO2 

-  —  (w/  -  OWI+/  +  “»+/)  (16) 

8 

The  wavenumber  k  is  averaged  in  the  y  direction,  and  is  a 
consequence  of  the  parabolic  approximation  (Radder,  1979), 
which  in  practical  terms  limits  the  obliquity  of  the  wave  with 
respect  to  the  x  axis  of  the  grid.  Details  of  the  derivation  are 
presented  in  Kaihatu  and  Kirby  (1995)  and  not  replicated  here. 
The  primary  advantage  of  the  model  is  that  the  shoaling  and 


Fig.  2.  Dissipation  D  (normalized  by  the  maximum)  as  a  function  of  relative  water  depth  kh.  Water  depth  h=  1  m.  Solid  line:  {=  10,  y  =  0.9,  dm= 0.02  m.  Dashed  line: 
C  =  100,  y=0.9,  dm= 0.2  m.  Dotted  line:  C=  10,  7=0.5,  dm= 0.02  m.  Dash-dot  line:  100,  7=0.5,  dm=0.2  m. 
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dispersion  are  determined  by  fiilly-dispersive  linear  theory,  thus 
overcoming  the  shallow  water  limitation  of  classical  Boussinesq 
type  models  (e.g.  Freilich  and  Guza  1984;  Liu  et  al.,  1985). 
Since  the  equation  is  parabolic,  numerical  solution  is  easily 
undertaken  via  the  Crank-Nicholson  method,  iterated  to  center 
the  nonlinear  terms  between  the  cross-shore  grid  points.  The 
process  of  parabolization  reduces  the  computational  effort,  but 
disallows  processes  which  run  backward  through  the  domain 
such  as  back  reflection.  Kaihatu  and  Kirby  (1995)  showed  that 
(14)  compared  well  to  laboratory  data.  The  parabolic  model  was 
extended  to  wide  angles  of  propagation  by  Kaihatu  (2001). 

Kaihatu  and  Kirby  (1995)  also  used  a  one-dimensional 
reduction  of  (14)  for  shoaling  waves,  and  extended  it  into  the 
surf  zone  by  the  addition  of  a  spectral  dissipation  term: 


4.  Model  results 

This  section  presents  a  number  of  consistency  and  model  data 
comparison  tests  performed  with  the  model  (18)  to  ascertain  the 
effect  of  the  mud  on  wave  transformation.  As  a  consistency 
check,  we  first  compare  the  model  to  laboratory  data  for  a  simple 
case  of  wave  propagation  over  a  flat,  muddy  bottom.  After 
confirming  the  performance  of  the  model,  we  then  use  it  to 
examine  the  physical  nature  of  nonlinear  wave  propagation  over 
viscous  muds.  In  particular,  we  analyze  the  effect  of  mud  on 
nonlinear  energy  transfer,  focusing  on  the  attenuation  of  high 
frequency  energy. 

4. 1.  Comparison  with  laboratory  data 
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where  a  is  the  dissipation  coefficient  due  to  surf  zone  dissipation; 
this  term  is  detailed  in  Kaihatu  and  Kirby  ( 1 995).  Comparisons  with 
data  yielded  excellent  results.  A  nonlinear  correction  term  was 
added  to  the  model  by  Kaihatu  (2001 ),  which  led  to  improved  wave 
shape  predictions.  We  do  not  employ  the  surf  zone  dissipation  in 
this  study,  but  reserve  it  for  analysis  in  a  future  study.  This  equation 
was  numerically  solved  with  a  fourth-order  Runge-Kutta  scheme. 

Similar  to  the  development  of  the  breaking  term  implementation 
in  Kirby  and  Dalrymplc  ( 1 986),  we  add  the  dissipation  mechanism 
of  Ng  (2000)  directly  to  the  equations  in  a  manner  consistent  with 
the  conservation  of  energy  flux.  Thus  (14)  becomes: 
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and  (17)  becomes: 
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Most  of  the  existing  field  or  laboratory  data  sets  on  wave 
propagation  over  muddy  sea  beds  (e.g.  Gadc  1958;  Tubman  and 
Suhayda,  1976;  Forristall  and  Recce,  1985;  Jiang  and  Mehta, 
1995)  focus  on  linear  processes  of  narrow-spectrum  long 
waves.  Their  approach  is  based  on  the  assumption  that  bottom 
interaction  is  significant  only  for  the  long  waves,  which  reach 
deep  enough  in  the  water  column  to  interact  directly  with  the 
bottom.  A  number  of  field  experiments  were  conducted 
recently,  dedicated  to  monitoring  nonlinear  aspects  of  wave 
evolution  (e.g.  Jaramillo  et  al.,  2006)  but  as  this  paper  was 
written  the  data  from  these  experiments  were  still  being 
processed.  A  detailed  analysis  of  the  performance  of  the 
model  in  field  data  tests  will  be  presented  elsewhere.  Flere,  we 
will  limit  scope  of  the  model-data  comparison  to  some  simple 
consistency  checks  against  laboratory  data  of  De  Wit  (1995), 
who  investigated  the  liquefaction  of  cohesive  sediments  in  a 
flume  under  varying  conditions  of  waves  and  currents. 

The  flume  used  to  perform  the  experiments  was  40  m  long, 
with  a  width  and  depth  of  0.8  m.  The  flume  was  fitted  with  a 
false  floor  of  depth  0.2  m  and  length  8  m,  allowing  an  overall 
maximum  depth  of  1  m.  In  this  area  dense  mud  was  emplaced. 
De  Wit  ( 1 995)  performed  tests  using  several  clays;  we  use  his 
results  from  Experiment  III,  in  which  China  clay  (kaolinite  is  a 
principal  mineral  component)  was  used.  The  clay  powder  was 
mixed  with  0.5%  solution  of  sodium  chloride  to  a  sediment 
concentration  of  275  kg/m3,  then  allowed  to  consolidate  for 
6  days.  This  resulted  in  a  suspension  with  a  density  of  about 
1300  kg/m3.  Its  viscosity  vm  =  2.7xl0  3  m2/s.  We  assumed 
v*,— 1.3  x  10” 6  m2/s. 

Two  other  tests  in  Experiment  III  were  carried  out  over  this 
material.  The  material  was  then  allowed  to  rest  for  one  day.  Test 
3  of  Experiment  III  was  then  conducted.  The  mud  parameters 
and  wave  conditions  used  are  shown  in  Table  1 .  This  experiment 
was  also  used  to  validate  a  more  comprehensive  mud  dissipation 
model  (Winterwerp  et  al.,  2007). 


Table  1 

Mud  parameters  and  wave  condition:  Test  3  of  Experiment  III  of  Dc  Wit  ( 1 995) 


pw  (kg/m3) 

Pm  (kg/m3)  dm  (m)  h  (m) 

H  (m)  v„  (m2/s) 

ns) 

1000 

1300  0.115  0.325 

0.045  2.6x1  O'5 

1.5 

where  Dn=D(conJclnJh)  is  as  in  (6). 
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Fig.  3.  Comparison  of  model  to  data  from  De  Wit  (1995),  Test  3  of  Experiment  III.  Solid  line:  waveheights  from  model.  Circles:  measurements. 


We  used  a  permanent  form  solution  of  (17)  without  the 
dissipation  terms.  This  solution  consisted  of  seven  harmonics 
phase-locked  to  the  primary  mode  with  amplitudes  spanning  six 
orders  of  magnitude;  details  on  this  solution  technique  are 
provided  by  Kaihatu  (2001).  The  amplitude  series  comprising 
this  waveform  were  input  into  the  model  and  run  over  the  8  m 
long  mud  section.  Waveheights  were  calculated  at  the  gage 


locations  by  reconstructing  the  waveform  and  determining  the 
distance  between  crest  and  trough.  The  resulting  comparison  is 
shown  in  Fig.  3;  agreement  is  good  and  on  a  par  with  that  shown 
by  Winterwerp  et  al.  (2007).  This  is  particularly  noteworthy  in 
that  the  depth  of  the  mud  layer  is  roughly  33%  of  the  water 
depth,  which  would  appear  to  violate  the  thin  layer  assumption 
inherent  in  the  formulation  of  Ng  (2000).  In  fact,  the  largest  of 


Fig.  4.  One-dimensional  cnoidal  waves  over  mud  patch:  {  =  10  case.  Solid  line:  7=0.25.  Dashed  line:  7=0.5.  Dash-dot  line:  7=0.75.  Dotted  line:  7  =  0.9. 
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the  scales  (2)  used  by  Ng  (2000)  is  kdm= 0.3,  which  is  not  very 
small.  However  the  fact  that  the  model  still  compared  well  to 
data  indicates  that  the  scaling  assumption  is  perhaps  not  that 
restrictive,  at  least  as  far  as  the  damping  is  concerned.  We  note 
here  that  the  sudden  increase  in  measured  waveheight  near  the 
downwave  end  of  the  mud  patch  was  likely  caused  by  wave 
reflection  from  the  back  of  the  mud  pit.  Due  to  the  parabolic 
nature  of  the  model,  this  effect  cannot  be  simulated.  Using 
linear  long  wave  theory  (Dean  and  Dalrymple,  1991)  the 
expected  reflection  is  estimated  to  be  10%  of  the  incident  wave 
energy.  This  slightly  overestimates  the  error  between  model  and 
measurement  beyond  jc=5  m,  possibly  due  to  the  use  of  linear 
long  wave  theory  for  its  estimation. 

4.2.  Permanent  form  shallow  water  waves  over  a  flat  bottom  — 
effects  of  mud 

Similar  to  the  previous  test,  our  next  test  involves  the 
transformation  of  waves  of  permanent  form  over  a  flat  bottom  in 
which  a  region  of  mud  has  been  placed.  However,  in  this  test  we 


will  investigate  the  damping  of  shallow  water  waves  across  long 
transformation  distances.  Our  domain  is  1000  m  long,  with  a 
constant  depth  of  1  m  and  a  grid  resolution  Ajc=0.025  m.  In 
actual  practice  this  is  far  finer  than  required  for  an  accurate 
solution,  but  is  used  to  sharply  resolve  the  free  surface  t]  in  a 
post-processing  step  (via  Eq.  (13))  once  the  complex  amplitudes 
An  are  calculated  by  the  model.  We  define  a  “mud  patch”  which 
extends  from  jc=300  m  to  jc=800  m.  We  note  here  that  the 
length  of  the  domain,  the  water  depth  and  the  extent  of  the  mud 
patch  were  chosen  to  bring  forth  particular  features  of  the  mud 
dissipation  effect  on  wave- wave  interaction,  as  will  be  seen. 
However,  these  values  are  comparable  to  many  areas  in  the 
world  with  long  shallow  shelves  and  muddy  bottoms.  For 
example,  the  shelf  near  Marsh  Island  in  Louisiana,  USA  is 
extremely  flat  and  less  than  2  m  in  depth  over  a  distance  of 
around  7  km  (Sheremet,  personal  communication).  As  before, 
reflection  of  the  wave  energy  is  not  considered  herein.  From 
linear  long  wave  theory  (Dean  and  Dalrymple,  1991)  the 
expected  reflection  due  to  the  mud  patch  is  estimated  to  be  less 
than  5%  of  the  incident  energy. 


Fig.  5.  One-dimensional  cnoidal  waves  over  mud  patch:  £  =  100  case.  Solid  line:  y=0.25.  Dashed  line:  y=0.5.  Dash-dot  line:  y=0.75.  Dotted  line:  7=0.9. 
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We  investigate  the  wave  behavior  with  two  different  sets  of 
parameters  to  describe  the  mud  layer.  The  first  uses  a  depth  of 
mud  of  0.02  m  and  a  mud  kinematic  viscosity  vm=  1.3  *  10”4, 
leading  to  £=  10  if  we  assume  the  kinematic  viscosity  of  water 
vw=1.3x  10~6.  The  dimensionless  mud  depth  d=dm/6m  =  1. 
The  second  case  uses  vm=  1.3  x  10~2,  making  C =  1 00.  In  order 
to  keep  the  dimensionless  mud  depth  d=  1  we  increase  the  mud 
layer  depth  dm=0.2  m.  We  will  simply  refer  to  them  as  the 
C=  10  case  and  the  £=100  case.  For  each  case  we  vary  the 
density  ratio  y= (0.25,0.5,0.75,0.9),  with  the  most  dampening 
occurring  for  y=0.9,  keeping  £  fixed. 

Our  shallow  water  incident  wave  is  a  permanent  form  solution 
of  (17),  with  a  period  T=  10  s  and  a  height  H=  0.1  m.  With  this 
wave  period,  we  find  that  kdm= 0.04,  which  is  sufficiently  small. 
The  solution  is  generated  via  the  technique  described  in  Kaihatu 
(2001 )  using  N=  15  frequency  components. 

Fig.  4  shows  free  surface  elevations  for  the  £  =  10  case.  As 
expected  the  greatest  amount  of  dampening  was  seen  when 
y=0.9.  However,  it  is  also  interesting  that  the  modification  to 
the  wavelength  by  the  presence  of  mud  did  not  cause  significant 
detuning  of  the  phase  locking  required  to  maintain  permanent 


form;  the  wave  retained  a  mostly-cnoidal  shape  within  and 
beyond  the  mud  patch  even  for  high  y  values. 

Fig.  5  shows  the  case  for  £  =  100.  The  damping  is  far  more 
dramatic  for  this  case  than  for  £  =  10;  for  the  y=0.9  case  the 
wave  almost  disappears  entirely.  To  the  extent  that  the  wave 
shape  is  discernible  at  high  y,  a  somewhat  more  sinusoidal 
shape  is  apparent. 

The  amplitudes  of  the  frequency  components  of  the 
wavefield  are  shown  in  Fig.  6  for  both  the  £=10  and  £=100 
cases,  for  all  values  of  y  used.  For  both  cases,  both  the  initial 
condition  at  jc=0  and  the  resulting  amplitudes  at  the  end  of  the 
mud  patch  (jc=800  m)  are  shown.  The  monotonic  decrease  in 
amplitude  with  increasing  frequency  is  apparent  for  the  £=10 
case  at  x=800  m,  and  is  also  reflected  in  the  retention  of  the 
cnoidal-like  shape  of  the  wavefield  (Fig.  4).  In  contrast,  the 
£=100  case  shows  a  distinct  non-monotonicity  with  frequency; 
the  fourth  harmonic  of  the  wave  (/=0.4  Hz)  displays  the  greatest 
amount  of  dissipation  (relative  to  the  initial  amplitude)  from  the 
mud  patch  in  the  range  /<  1  Hz  (kh  ~tt).  At  /<0.4  Hz,  kh 
~0.89,  which  is  close  to  the  point  of  maximum  dissipation 
displayed  in  Fig.  2. 


Fig.  6.  One-dimensional  cnoidal  waves  over  mud  patch:  amplitudes  of  wavefield  at  jc=0  (dash-x  line)  and  at  downwave  end  of  mud  patch  (jc=800  m).  Solid  line: 
y=0.25.  Dashed  line:  y=0.5.  Dash-dot  line:  y=0.75.  Dotted  line:  y=0.9.  Top:  £=10.  Bottom:  £=100. 
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Interestingly,  it  is  also  evident  that  there  is  a  significant  loss 
of  amplitude  at  the  high  frequencies  over  the  mud  patch.  The 
relative  depth  kh  is  around  9  at  /=  1.5  Hz,  far  beyond  the 
kinematical  deep  water  limit.  It  is  thus  likely  that  subharmonic 
interaction  with  the  lower,  actively-damped  frequencies  is 
causing  this  reduction  in  amplitude,  as  shown  by  Shcremet  et  al. 
(2005)  for  the  case  of  a  wave  group.  To  investigate  this  we  re¬ 
run  the  model  for  the  case  of  7=0.9,  using  both  £=10  and 
£=100.  The  first  run  is  linear,  with  the  nonlinearity  in  (17) 
deactivated.  The  second  run  reduces  the  nonlinear  interactions 
by  deactivating  the  subharmonic  interactions  (the  first  summa¬ 
tion  on  the  right  side  of  Eq.  (17)).  If  the  subharmonic  inter¬ 
actions  are  in  fact  causing  the  depletion  of  energy  in  frequencies 
beyond  the  deep  water  limit,  we  would  expect  to  see  similar 
amplitudes  between  the  linear  and  reduced  nonlinear  models. 

Fig.  7  compares  the  resulting  evolution  characteristics 
between  the  full  model  (17),  the  linearized  model,  and  the 
reduced  nonlinear  model  (with  subharmonic  evolution  deacti¬ 
vated),  for  y=0.9  and  both  £=10  and  £=100.  There  is  a  close 
correspondence  between  the  linear  and  the  reduced  nonlinear 
model  at  the  downwave  side  of  the  mud  patch,  indicating  that 
the  (deactivated)  subharmonic  interactions  are  responsible  for 


the  reduction  in  amplitudes  at  high  frequencies,  as  seen  in  the 
results  of  the  full  model.  The  high  frequency  amplitudes  of  the 
reduced  nonlinear  model  are  elevated  slightly  over  those  of  the 
linear  model  due  to  superharmonic  interactions  moving  energy 
from  the  low  frequencies. 

5.  Two  dimensional  wave  propagation  over  mud 

Our  next  set  of  simulations  involves  the  propagation  of  a  two 
dimensional  wave  field  over  a  small  patch  of  mud,  with 
emphasis  placed  on  the  refraction  and  diffraction  patterns 
apparent  over,  and  in  the  lee  of,  the  patch. 

We  use  a  two-dimensional  domain  that  is  2000  m  in  the 
propagation  direction  by  600  m  in  the  longshore  direction,  and 
use  closed  lateral  boundaries.  The  extent  in  the  propagation 
direction  should  allow  enough  distance  to  display  significant 
diffraction  effects  in  the  lee  of  the  mud  patch.  We  use  the 
T=  10  s  wave  used  previously  for  the  one-dimensional  case  as 
our  initial  condition.  The  bathymetry  is  a  constant  depth 
(h=  1  m)  as  before.  Into  this  domain  we  add  a  region  of  mud 
(200<jcmMrf<500  m,  150  m<ymud <450  m)  with  dm- 0.2  m, 
7  =  0.9  and  £=100.  The  grid  resolution  used  is  Ax =0.3  m  and 


Fig.  7.  One-dimensional  cnoidal  waves  over  mud  patch  with  y=0.9:  amplitudes  of  wavefield  at  jc=0  (dash  -x  line)  and  at  downwave  end  of  mud  patch  (*=800  m). 
Solid  line:  Eq.  (1 1).  Dashed  line:  Linear  model.  Dash-dot  line:  Eq.  (1 1)  with  only  superharmonic  interactions.  Top:  £=10.  Bottom:  C  =  100. 
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Fig.  9.  Cnoidal  waves  propagating  over  two-dimensional  mud  patch:  overhead  view  of  free  surface  elevation  (crests  are  in  light  streaks).  Mud  patch  boundary  shown  in 
dashed  lines,  and  is  reduced  in  longshore  extent.  Mud  parameters:  dm= 0.2  m,  y=0.9,  {=  100. 
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A>>=0.6  m;  again  this  is  far  finer  than  required  for  accurate 
solution  but  serves  to  help  sharpen  the  depiction  of  the  free  surface 
from  the  complex  amplitudes.  (We  note  here  that  our  resolution 
in  the  propagation  direction  Ax  is  decreased  relative  to  that  of 
the  one-dimensional  case  due  to  the  extra  memory  requirements 
incurred  by  the  expansion  to  two  horizontal  dimensions.) 

The  result  is  shown  in  Fig.  8.  There  is  a  very  strong  attenuation 
of  the  wave  elevation  through  the  mud  patch  and  a  very  long 
shadow  zone  in  the  lee  of  the  patch.  Additionally,  the  wave  crests 
inside  the  mud  patch  appear  rounder  and  less  sharp  than  those 
outside  the  patch,  likely  due  to  the  preferential  dampening  of  the 
higher  frequency  components  of  the  wave.  As  the  waves 
propagate  past  the  mud  patch,  wave  diffraction  moves  energy 
into  the  shadow  region.  Due  to  the  effect  of  mud  on  the 
wavelength  inside  the  mud  region,  the  wave  crests  in  the  lee  of  the 
region  become  discontinuous  with  those  outside  it;  this  effect  is 
also  seen  in  wave  propagation  over  two-dimensional  varying 
bathymetry  (e.g.,  Berkhoff  ct  al.,  1982).  Eventually,  further 
downwave  of  the  mud  region,  the  bisected  wave  crests  intersect, 
resulting  in  a  hexagonal  pattern  similar  to  the  intersecting  cnoidal 
waves  generated  in  the  laboratory  by  Hammack  et  al.  (1989). 

In  the  classic  case  of  wave  diffraction  around  a  detached 
breakwater,  the  extent  of  the  shadow  zone  is  strongly  dependent 
on  the  width  of  the  breakwater  (Penney  and  Price,  1952).  We 
anticipate  that  a  similar  dependence  is  evident  on  the  longshore 
dimension  of  the  mud  patch.  We  decrease  the  longshore  extent 
of  the  mud  patch  considerably  (250  rn<j>mud<350  m)  and 
rerun  the  simulation.  The  result  is  shown  in  Fig.  9.  Similar  to  the 
problem  of  a  detached  breakwater,  the  longshore  extent  of  the 
shadow  zone  in  the  lee  of  the  mud  patch  is  closely  related  to  that 
of  the  mud  patch.  In  contrast  to  Fig.  8,  here  we  see  a  much 
narrower  shadow  zone.  Additionally,  well  into  the  lee,  the 
intersection  of  the  crests  generate  a  wave  traveling  at  the  mean 
direction,  a  nonlinear  effect  similar  to  the  hexagonal  patterns 
seen  in  Fig.  8.  From  the  point  of  view  of  using  remote  sensing  to 
seed  inverse  models  to  infer  the  nature  of  the  bottom,  this 
sensitivity  to  patch  size  is  encouraging  provided  there  is 
sufficient  damping  of  the  waves  by  the  bottom  mud. 

6.  Summary 

In  this  study  we  have  outlined  the  development  of  a  nonlinear 
wave  model  adapted  with  a  mechanism  for  dissipation  by  a  thin 
layer  of  viscous  mud.  The  viscous  damping  mechanism  used 
was  formulated  by  Ng  (2000)  as  a  boundary-layer  reduction  of 
the  model  of  Dalrymplc  and  Liu  ( 1 978)  for  damping  by  a  lower 
layer  of  viscous  mud  of  arbitrary  depth. 

The  dispersion  relation  (5)  resulting  from  the  theory  of  Ng 
(2000)  is  complex;  the  rate  of  viscous  damping  is  contained  in 
the  imaginary  part  of  the  dispersion  relation  (6),  while  the  effect 
of  the  mud  on  wave  kinematics  is  contained  in  the  real  part  (12). 
Ng  (2000)  detailed  the  damping  effect  of  the  theory  but  did  not 
investigate  the  effect  on  the  wavelength.  We  showed  that  the 
theory  of  Ng  (2000)  does  decrease  the  wavenumber  as  d ,  the 
normalized  mud  depth,  increases;  however,  the  rate  of  this 
decrease  is  less  than  that  shown  by  Dalrymplc  and  Liu  (1978). 
Thus,  the  viscous  mud  layer  will  not  alter  the  detuning  of  the 


nonlinear  interactions  (dependent  on  the  wavenumbers  of  the 
interacting  triads)  relative  to  that  with  no  mud  present.  The 
dependence  of  the  damping  rate  (6)  on  the  relative  water  depth 
kh  is  shown  in  Fig.  2;  there  is  an  apparent  maximum  damping 
near  kh  ~0.88,  though  the  location  of  this  maximum  will 
change  with  different  h  and  d. 

The  parabolic  nonlinear  model  of  Kaihatu  and  Kirby  (1995) 
is  used  as  the  primary  wave  model  for  weakly  nonlinear, 
fully  dispersive  waves,  with  the  implementation  of  the  viscous 
damping  mechanism  of  Ng  (2000)  performed  in  a  manner 
similar  to  the  surf  zone  dissipation  model  of  Kirby  and 
Dalrymple  (1986).  The  model  was  then  tested  against  the 
experimental  data  of  Dc  Wit  ( 1 995)  and  shown  to  compare  well, 
despite  the  thin  layer  assumption  behind  the  dissipation 
mechanism. 

Further  testing  was  then  performed  in  both  one  and  two 
spatial  dimensions,  consisting  of  a  permanent  form  cnoidal 
wave  propagating  over  mud  patches.  Two  general  cases  for  mud 
viscosity  were  detailed.  The  first  used  a  mud  viscosity 
vm=1.3xl(T4  m2/s,  or  a  viscosity  ratio  (=10.  The  second 
used  vm=  1.3  x  10“  2  m2/s,  or  C  =  100.  For  each  case  the  density 
ratio  y  = (0.25,0.5,0.75,0.9).  The  mud  depth  dm  was  selected  for 
each  viscosity  case  such  that  the  dimensionless  ratio  d=dj 
<5m=l,  where  <5m  is  the  Stokes  boundary  layer  thickness  in  the 
mud  layer. 

For  the  one-dimensional  scenario,  the  cnoidal  waveform 
retained  its  general  shape  for  the  (=10  case,  while  for  the 
highly-damped  (=100  case  the  waveform  (to  the  extent  that 
dissipation  did  not  obscure  it)  appeared  to  be  more  sinusoidal. 
This  indicates  that  there  is  some  preferential  dampening  at  the 
higher  frequencies  of  the  wavefield,  in  line  with  the  kh- 
dependent  dissipation  shown  in  Fig.  2. 

To  elucidate  this  further,  the  amplitudes  of  the  wavefield 
were  plotted  for  the  different  y  values  for  each  viscosity  case 
(Fig.  6),  both  at  the  initial  condition  and  at  the  downwave  side 
of  the  mud  patch.  While  the  amplitudes  for  the  (=10  case 
monotonically  decreased  with  frequency  at  the  end  of  the  mud 
patch,  the  (=  100  case  did  not,  with  the  fourth  harmonic  of  the 
wavefield  exhibiting  the  greatest  dissipation  relative  to  the 
initial  condition.  For  this  harmonic  A:/i  —  0.89,  which  is  close  to 
the  location  of  the  dissipation  maximum  shown  in  Fig.  2. 

It  was  also  apparent  that  the  frequencies  beyond  the  deep 
water  limit  displayed  some  degree  of  attenuation  over  the  mud 
patch  for  both  viscosity  cases,  perhaps  due  to  subharmonic 
interactions  between  them  and  lower,  actively-dissipated 
frequencies.  The  one-dimensional  model  was  run  again  for 
both  (=10  and  (=100  (y  =0.9).  One  run  was  performed  with 
nonlinearity  deactivated,  and  another  with  subharmonic  inter¬ 
actions  (first  summation  on  right  side  of  Eq.  (17))  deactivated.  It 
was  shown  (Fig.  7)  that  the  amplitudes  for  both  the  linear  and 
reduced  nonlinear  cases  matched  very  closely  and  exhibited  no 
reduction  from  the  mud  patch  in  the  high  frequencies.  This 
indicates  that  the  (deactivated)  subharmonic  interactions 
are  responsible  for  the  dissipation  of  the  high  frequencies  of 
the  wavefield,  even  if  beyond  the  deep  water  limit.  This  is  in 
line  with  the  finding  of  Sheremet  et  al.  (2005)  for  the  case  of 
wave  groups. 
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The  two-dimensional  model  was  then  applied  to  investigate 
wave  propagation  over  a  mud  patch  of  limited  longshore  extent. 
Fig.  8  shows  a  result  for  a  mud  patch  whose  longshore 
dimension  is  half  that  of  the  domain.  Very  strong  damping 
occurs  inside  the  patch,  leading  to  significant  diffraction  on  the 
lee  side.  Further  down  wave,  the  discontinuity  in  the  wave  crest 
caused  by  the  mud  patch  begins  to  resolve  into  a  hexagonal 
wave  pattern,  as  the  two  halves  of  the  bisected  wave  crest  rejoin. 
Limiting  the  longshore  extent  of  the  mud  patch  (Fig.  9)  leads  to 
a  very  different  pattern  on  the  downwave  side.  Since  less  of  the 
wave  crest  experiences  damping,  there  is  less  diffraction. 

Future  work  with  the  model  will  involve  surf  zone 
dissipation  and  the  resultant  damping  compared  to  that  caused 
by  mud.  Of  particular  interest  would  be  the  alteration  of  the 
velocity  moments  required  for  nearshore  sediment  transport 
calculations.  We  also  anticipate  the  opportunity  to  perform 
comparisons  with  field  data. 
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